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We consider the problem of determining whether genuine multipartite entanglement was produced 
in an experiment, without relying on a characterization of the systems observed or of the measure- 
ments performed. We present an n-partite inequality that is satisfied by all correlations produced 
by measurements on biseparable quantum states, but which can be violated by n-partite entangled 
states, such as GHZ states. In contrast to traditional entanglement witnesses, the violation of this 
inequality implies that the state is not biseparable independently of the Hilbert space dimension 
and of the measured operators. Violation of this inequality does not imply, however, genuine mul- 
tipartite non-locality. We show more generically how the problem of identifying genuine tripartite 
entanglement in a device-independent way can be addressed through semidefinite programming. 



The generation of multipartite entanglement is a cen- 
tral objective in experimental quantum physics. For in- 
stance, entangled states of fourteen ions and six photons 
have recently been produced [H [5] . In any such experi- 
ment, a typical question arises: How can we be sure that 
genuine n-partite entanglement was present? A state is 
said to be genuinely n-partite entangled if it is not bisep- 
arable, that is, if it cannot be prepared by mixing states 
that are separable with respect to some partition. Con- 
sider for instance the tripartite case: a state pbs is said 
to be biseparable if it admits a decomposition 

^ bs = 12pab®Pc + Y1pac®Pb+ S ^Pbc®Pai (1) 

where the weight of each individual state in the mix- 
ture has been included in its normalization; a state that 
cannot be written as above is genuinely tripartite entan- 
gled. Determining whether genuine n-partite entangle- 
ment was produced in an experiment represents a diffi- 
cult problem that has attracted much attention recently 
(see e.g. [31 2]). The usual approach consists of measur- 
ing a witness of genuine multipartite entanglement, or 
of doing the full tomography of the state followed by a 
direct analysis of the reconstructed density matrix. 

Such approaches, however, not only rely on the ob- 
served statistics to conclude about the presence of en- 
tanglement, but also require a detailed characterization 
of the systems observed and of the measurements per- 
formed. Consider for instance the following witness of 
genuine tripartite entanglement: 

M = X X X 2 X 3 - XxY 2 Y 3 - Y X X 2 Y Z - Y 1 Y 2 X 3 , (2) 

where Xj — a x and Yj — a y are the Pauli spin observ- 
ables in the x and y direction for particle j. For any 
biseparable three-qubit state (M) = tr (Mp) < 2 [5] • 
Thus if we measure three spin-| particles in the x and 
y direction and find an average value (M) > 2, we can 
conclude that the state exhibits genuine tripartite entan- 
glement. 

Suppose, however, that the measurement Ys carries a 
slight (possibly unnoticed) bias towards the x direction. 



That is, instead of measuring Y 3 = a y , we actually mea- 
sure Y 3 — cos9a y +sm9a x . Then it is not difficult to see, 
all other measurements being ideal, that the biseparable 
state \f) = \ (|00) + e-^\l\)) AB ® (|0) + |l)) c , where 

<f> = arctan(sin#), yields (M) = 2 \J\ + sin 2 9 which is 
strictly larger than 2 for any 9^0. Thus, unless we 
measure all particles exactly along the x and y directions, 
we can no longer conclude that observing (M) > 2 im- 
plies genuine tripartite entanglement. Importantly, this 
is not a unique feature of the above witness, but rather 
all conventional witnesses are, to some extent, suscepti- 
ble to such systematic errors that are seldom taken into 
account. 

Furthermore, tomography and usual entanglement wit- 
nesses typically assume that the dimension of the Hilbert 
space is known. For instance, in a typical experiment 
demonstrating, say, entanglement between four ions, we 
usually view each ion as a two-level system. But an ion 
is a complex object with many degrees of freedom (posi- 
tion, vibrational modes, internal energy levels, etc). How 
do we know, given the inevitable imperfections of the ex- 
periment, that it is justified to treat the relevant Hilbert 
space of each ion as two-dimensional and how does this 
simplification affects our conclusions about the entangle- 
ment present in the system [5]? Even if it is justified to 
view each ion as a qubit, is entanglement between four 
systems really necessary to reproduce the measurement 
data, or could they be reproduced with fewer entangled 
systems if qutrits were manipulated instead? 

These remarks motivate the introduction of entangle- 
ment witnesses that are able to guarantee that a quan- 
tum system exhibits (n-partite) entanglement, without 
relying on the types of measurements performed, the 
precision involved in their implementation, or on as- 
sumptions about the relevant Hilbert space dimension. 
We call such witnesses, device-independent entanglement 
witnesses (DIEW). This type of approach was already 
considered in Refs. [7J[5]. Note that other solutions to 
the above problems are possible, such as entanglement 
witnesses tolerating a certain misalignement in the mea- 
surement apparatuses [5] or the characterization of real- 
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istic measurement apparatuses through squashing maps 
[6]. These types of more specific approaches, however, 
still require some partial characterization of the system 
and measurement apparatuses, which is not necessary 
when using DIEWs. 

Any DIEW is a Bell inequality (i.e. a witness of non- 
locality). Indeed, i) the violation of a Bell inequality 
implies the presence of entanglement, and ii) any mea- 
surement data that does not violate any Bell inequality 
can be reproduced using quantum states that are fully 
separable [TU]. The violation of a Bell inequality is thus 
a necessary and sufficient condition for the detection of 
entanglement in a device-independent (DI) setting. This 
observation is the main insight behind DI quantum cryp- 
tography [TUE2], where the presence of entanglement is 
the basis of security. 

The relation between DIEW for genuine n-partite en- 
tanglement and witnesses of multipartite nonlocality is 
more subtle. While there exist Bell inequalities that de- 
tect genuine n-partite nonlocality [T3J [TU [TU , not every 
DIEW for n-partite entanglement corresponds to such 
a Bell inequality. Consider for instance, the expression 
([2| . If no assumptions are made on the type of systems 
observed and measurements performed, the inequality 
(M) < 2 corresponds to Mermin's Bell-type inequality 
[17] . i.e., a value (M) > 2 necessarily reveal non- locality, 
hence entanglement. Moreover, a value (M) > 2-\/2 guar- 
antees genuine tripartite entanglement [7J [14] . The Mer- 
min expression |2]) can thus be used as a tripartite DIEW. 
Yet, it cannot be used as a Bell inequality for genuine tri- 
partite non-locality, since a simple model involving com- 
munication between two parties only already achieves the 
algebraic maximum (M) = 4 [T"4l . 

The objectives of this paper are to formalize the con- 
cept of DIEW for genuine multipartite entanglement and 
initiate a systematic study that goes beyond the early 
examples given in [7J [5]. Following this line, we start 
by introducing the notion of quantum biseparable corre- 
lations. We then present a simple DIEW for n-partite 
entanglement which is stronger for GHZ (Greenberger- 
Horne-Zeilinger) states than all the inequalities intro- 
duced in [7J [8]. In the case n = 3, we also provide a 
general method for determining whether given correla- 
tions reveal genuine tripartite entanglement and apply 
it to GHZ and W states. Apart from yielding practical 
criteria for the characterization of entanglement in a mul- 
tipartite setting, our results also clarify the relation be- 
tween device-independent multipartite entanglement and 
mulitpartite nonlocality. 

1. Biseparable quantum correlations. For sim- 
plicity of exposition, let us consider an arbitrary tripar- 
tite system (the following discussion easily generalizes to 
the n-party case). To characterize in a DI way its en- 
tanglement properties, we consider a Bell-type experi- 
ment: on each subsystem, one of m possible measure- 
ments is performed, yielding one of d possible outcomes. 
We adopt a black-box description of the experiment and 
represent the measurements on each of the three sub- 



systems by classical labels x,y,z E {1, ...,m} (corre- 
sponding, e.g., to the values of macroscopic knobs on the 
measurement apparatuses) and denote the corresponding 
classical outcomes a,b,c € {1, . . . , d}. The correlations 
obtained in the experiment are characterized by the joint 
probabilities P(abc\xyz) of finding the triple of outcomes 
a, 6, c given the measurement settings x, y, z. 

We say that the correlations P(abc\xyz) are bisepa- 
rable quantum correlations if they can be reproduced 
through local measurements on a biseparable state pbs- 
That is, if there exist a biseparable quantum state ([I]) 
in some Hilbert space H, measurement operators M a \ x , 
Muy, and M c \ z (which without loss of generality we can 
take to be projections satisfying M a \ x M a ,\ x = S a ^M a \ x 
and J2 a M a \ x = 1), such that 



P(abc\xyz) = tr [M a \ x ® M b \ y <g> M r \ z p hs 



(3) 



If given quantum correlations P(abc\xyz) are not bisep- 
arable, they necessarily arise from measurements on a 
genuinely tripartite entangled state, and this conclusion 
is independent of any assumptions on the type of mea- 
surements performed or on the Hilbert space dimension. 

Equivalently, biseparable quantum correlations can be 
defined as those that can be written in the form 



P(abc\xyz) =Y,PQ(ab\xy)P%(c\z) 

k 

+ £ P*{ac\xz)P k Q {b\y) + P k Q {bc\yz)P*{a\x), 



(4) 



where Pg(a6|xy) and Pq(c\z) correspond, respectively, to 
arbitrary two-party and one-party quantum correlations, 
i.e., they are of the form Pg(ab\xy) = te[M$ x <8)Mk p k AB \ 

and Pq{c\z) = tr[M*j z pp] for some unormalized quan- 
tum states p\ B i Pq and measurement operators M*< x , 
M^ y , M^ z (and similarly for the other terms in (pit)). 
Note that here the measurement operators for different 
bi-partitions do not need to be the same (though this can 
always be achieved as shown in Appendix Dj. That any 
correlations of the form ^ is of the form |4 ) is immediate 
using the definition ([I]) of biseparable states. Conversely, 
it is easy to see that any correlations of the form Q is 
of the form see Appendix [A] 

Let Q3 denote the set of tripartite quantum corre- 
lations and Q2/1 C Q3 the set of biseparable quan- 
tum correlations. From |4j), it is clear that Q2/1 is 
convex and that its extremal points are of the form 
Pq xt (ab\xy)PQ Xt (c\z) where Pq :t (ab\xy) is an extremal 
point of the set Q2 of bipartite quantum correlations 
and Pg xt (c\z) an extremal point of the set Q\ of single- 
party correlations (the extreme points of Q\ are ac- 
tually classical, deterministic points). Since the set 
Q2/1 is convex, it can be entirely characterized by lin- 
ear inequalities. Those linear inequalities separating 
Q2/1 from Q 3 correspond to DIEWs for genuine tri- 
partite entanglement. Since Q2/1 has an infinite num- 
ber of extremal points, there exist an infinite number 
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FIG. 1. A particular slice of the space of tripartite correlations 
with 3 settings and 2 outcomes representing schematically the 
sets of general quantum correlations (Q3), Svetlichny corre- 
lations (S'2/1) and biseparable quantum correlations (Q2/i)- 
The point 1 corresponds to random correlations and P to the 
GHZ correlations maximally violating the DIEW j5j), which 
is represented by the straight line I3; note that a DIEW can 
be violated by Svetlichny-local correlations. (The Svetlichny 
polytope S'2/1 can be determined exactly using linear pro- 
gramming, while Q3 and Q2/1 can be approximated efficiently 
using SDP techniques, see main text.) 



of such inequalities. Note that the set of local corre- 
lations P(abc\xyz) = J^k P k {a\x)P k (b\y)P k (c\z) is con- 
tained in Q2/x- This implies that any DIEW for gen- 
uine tripartite entanglement is a Bell inequality (though 
not necessarily a tight Bell inequality). Note also that 
the decomposition Q corresponds to a Svetlichny-type 
decomposition |13j where all bipartite factors are re- 
stricted to be quantum, whereas less restrictive con- 
straints (or even none in Svetlichny original definition 
P(abc\xyz) — ^ k P k {ab\xy)P k {c\z) + . . . ) are imposed 
on these bipartite terms in the definitions of multipar- 
tite non- locality [T3] ■ It follows that the set of genuinely 
bipartite non-local correlations is larger than the set of 
biseparable quantum correlations as illustrated in Fig. [l] 
Thus while any Bell inequality detecting genuine tripar- 
tite non-locality is a DIEW for genuine tripartite entan- 
glement, the converse is not necessarily true. All these 
observations extend to the n-party case. 

2. A DIEW for n-partite entanglement. We 
now present a DIEW for n parties, where each party i 
performs a measurement x% € {1,2,3} and obtains an 
outcome a% € { — 1,1}. We denote E{x) the correlator 
associated to the measurement settings x = (x\ . . . , x n ), 
i.e. the expectation value E{x) = J2a P{o>\x) n"=i a «' 
where a = (ai . . . , a n ) denotes an n-tuple of outcomes. 
Let E k — J2 S $(527=1 Xi ~ k)E(x) be the sum of corre- 
lators E(x) for which the measurement settings Xi of the 
n parties sum up to k. Let fk be a function such that 
fk+3 = — fk and taking successively the values [1,1,0] on 
the integers k = 0,1,2. Then the inequality 

3n 

I n = Y.h-nE k <2xT- z ' 2 (5) 

k—n 



is satisfied by all biseparable quantum correlations, and 
is thus a DIEW for genuine n-partite entanglement. The 
proof of this statement is based on the decomposition Q 
and is given in Appendix [B] The Svetlichny bound asso- 
ciated to the expression /„, on the other hand, is easily 
found to be 4 x 3™~ 2 > 2 x 3™~ 3 / 2 (see Appendix[l|) [24]. 

We now illustrate how this DIEW can be used to de- 
tect genuine multipartite entanglement. For this, let us 
consider the noisy GHZ state p = V\GHZ) nn {GHZ\ + 
(1 — V)l/2 n characterized by the visibility V. Carrying 
out the measurements cos((^j^ — g^K) a x +sin((^pi — 
^-)tt) <J y on all parties, we obtain I n = "& n ~ x l 2 V , which 
violates (jS) provided that V > 2/3. The DIEW ^ can 
thus detect in a DI way genuine n-partite entanglement 
in a noisy GHZ state for visibilities as low as V = 2/3. 
This significantly improves over the threshold visibility 
V = l/y/2 required to violate the DIEW based on the 
Mermin expression ^ or the different inequalities intro- 
duced in [8]. Note that in a DI setting it is not possible 
to detect the genuine multipartite entanglement in tri- 
partite GHZ states below V — 1/2 using projective mea- 
surements. Indeed, we have shown in this case that there 
exists a biseparable model reproducing all GHZ correla- 
tions (see Appendix [C|) . 

In the case n = 3, the DIEW ^ takes the form I3 = 

^3+^3-^3-^3+^3 < 6 ^3- lt therefore involves only 
18 expectation values, compared to 27 for a full tomogra- 
phy of a three-qubit system. Let us stress, however, that 
contrary to usual entanglement witnesses I3 is not re- 
stricted to two-dimensional Hilbert spaces, even though 
it uses observables with binary outcomes. For instance, if 
all parties perform the measurements 2\<fi(xi)\4>{xi) \ — 1 
with \<j>{xi)) = ^= (|0) + e i ( 6 ^-7W 18 |l}) on the three- 

qutrit state \ (|000) + | 111) + |222», then I 3 = 6^ + 
8/3, showing that the state is genuinely tripartite- 
entangled. 

3. General characterization of biseparable 
quantum correlations in the case n = 3. Though 
the DIEW ([5]) seems particularly well adapted to GHZ 
states, we cannot expect a single nor a finite set of DIEW 
to completely characterize the biseparable region, as il- 
lustrated in Fig. 1. It is thus desirable to derive a gen- 
eral method to decide whether arbitrary correlations are 
biseparable. Here we show how the semi-definite pro- 
gramming (SDP) techniques introduced in (THJ, [TH] can 
be used to certify that the correlations observed in an 
experiment are genuinely tripartite entangled. 

Our approach is based on the observation that the ten- 
sor product separation pab ® Pc at the level of states in 
the definition ([I]) of biseparable states can be replaced by 
a commutation relation at the level of operators. Specif- 
ically, let s = {AB/C,AC/B,BC/A} denotes the three 
possible partitions of the parties into two groups. Then, 
P{abc\xyz) are biseparable quantum correlations if and 
only if there exist three arbitrary (not necessarily bisepa- 
rable) states p s and three sets of measurement operators 
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State Vmin with two settings Vmi n with three settings 

GHZ 0.7071 ^1/^2 0.6667 ~ 2/3 

W 0.7500 ~ 3/4 0.7158 

TABLE I. Summary of numerical investigations. 



{^v M %> M c y such that 



P(abc\xyz) 



MK\ x ®M° lv ®M s clzP s 



(6) 



where measurement operators corresponding to an iso- 
lated party commute, i.e., [M^f ,M^J°] — 0, and 
similarly for the other partitions. The equivalence 
between ^ and ([6| is established in Appendix [d] 
The problem of determining whether given correlations 
P{abc\xyz) are biseparable thus amounts to finding a 
set of operators satisfying a finite number of algebraic 
relations (the projection defining relations of the type 
M s , M s ,, = 5 a a ,M s , , V M s , = 1 and the commuta- 

a\x a'\x a > a a\x' t-^ia a\x 

tion relations mentioned above) such that (|6| holds. Such 
a problem is a typical instance of the SDP approach intro- 
duced in (THUS] (see details in Appendix [d|). Specifically, 
it follows from the results of [TH] that one can define an 
infinite hierarchy of criteria that are necessarily satisfied 
by any correlations of the form (|6| and which can be 
tested using SDP. If given correlations do not satisfy one 
of these criteria, we can conclude that they reveal gen- 
uinely tripartite entanglement. Further, it is possible in 
this case to derive an associated DIEW from the solution 
of the dual SDP. Modulo a technical assumption, it can 
be shown that the hierarchy of SDP criteria is complete, 
that is, if given correlations are not biseparable this will 
necessarily show-up at some finite step in the hierarchy. 

4. Application to GHZ and W states. Using fi- 
nite levels of this hierarchy and optimizing over the pos- 
sible measurements, we investigated the minimal visibili- 
ties above which the GHZ state \GHZ) and the W state 
\W) = (|001) + |010) + |100))/-y/3 exhibit correlations that 
are not biseparable (and thus reveal genuine tripartite 
entanglement) in the case of two and three measurement 
settings per party. Our results are summarized in Table]!] 



For GHZ states, the reported visibility Knin = 2/3 for 
three measurements per party correspond to the thresh- 
old required to violate the DIEW pH) , suggesting that this 
DIEW is optimal in this case. In the case of two measure- 
ments per party, we could not lower the visibility below 
the threshold V m i n = l/v^2, which corresponds to the 
visibility required to violate the DIEW based on Mer- 
min expression ^ and the DIEWs introduced in [3|S]. 
Note, however, that for V > 1/ \2 the GHZ state violates 
Svetlichny's inequality [13] and thus exhibits genuine tri- 
partite nonlocality. Thus for GHZ state the DIEWs in- 
troduced in [8] do not improve over what can already 
be concluded using the standard notion of tripartite non- 
locality. On the other hand, our numerical explorations 
suggest that the visibilities V m i n = 2/3 for GHZ states 
with three measurements and V m [ n = 3/4 for W states 
with two measurements cannot be attained using the no- 
tion of genuine tripartite non-locality, illustrating the in- 
terest of the weaker notion of DIEW. 

Dicussion. To conclude, we comment on some possi- 
ble directions for future research. First of all, note that 
by identifying the measurement settings "Xi" with Xi — 1 
and "Yi" with Xi = 2, the two-setting DIEW based on 
Mermin expression ([2| can be written as E\ — E\ < 2\/2, 
which is of the same general form as the three-setting 
DIEW ([5|. This suggests that the DIEWs based on 
([2]) and ([5]) actually form part of a larger family of m- 
settings DIEWs. This question deserves further inves- 
tigation. A second problem is to derive simple DIEWs 
that are adapted to W states and that can in particular 
reproduce the threshold visibilities obtained in Table 1. 
Finally, we have shown a practical method to charac- 
terize three-partite biseparable correlations using SDP. 
It would be interesting to understand how to generalize 
these results to the n-partite case. A possibility would 
be to combine the approach of [TH with the symmet- 
ric extensions introduced in 20 . This question will be 
investigated elsewhere. 
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Appendix A: Biseparable quantum correlations 

Here, we show that the definitions ^ and (4]) of bisep- 
arable correlations are equivalent. For the sake of gener- 
ality, we prove this equivalence in the n-partite scenario. 
Let t denote a subset of {1, . . . , n} containing from one 
to n — 1 elements, and let t' = {1, . . . , n} \ t denote the 
complementary subset. The pair (t, t') thus represent a 
partition of the n parties into two non-empty groups. Let 
Pt <8> Pt' be an (unormalized) n-partite state that is prod- 
uct across the partition (i, t'). An n-partite state pb s is 
then biseparable if it can be written as a mixture 



Pbs 



(Al) 



k t 



A state that cannot be written in the above form is gen- 
uinely n-partite entangled. 

Let P(a\ . . . a n \x\ . . . x n ) — P{a\x) be an n-partite 
probability distribution characterizing a Bell experiment 
with measurement settings x = (x\ . . . x n ) and measure- 
ment outcomes a = (ai...a n ). We say that the cor- 
relations described by P(a\x) are biseparable quantum 
correlations iff 



P{a\x)=tr [M ailxi 



M, 



a n \x n Pbs\ 



(A2) 



for some biseparable state pi, s and measurement oper- 



ators M. 



CLi \Xi 



Equivalently, as we will see, biseparable 



quantum correlations can be defined as those that can 
be written in the form 



P{a\x) = Y,Y, P Q^) P Q^'\ 2 v 



(A3) 



fet 



where we write a t and Xt for the outcome and measure- 
ment settings of the parties belonging to the subset t 



and where Pg(at\xt) = tr 



are arbitrary 



(unnormalized) quantum correlations for the parties in t. 
That any correlations of the form (A2) is of the form 



( A3 ) is immediate using the definition (All of biseparable 



states. Conversely, (A3) can be rewritten in the form 



(A4) 



gives a representation that is almost of the 
( A2 1 , except that a different set of measure- 



This 
form 

ment operators corresponds to each term in the de- 
composition (Al) of the biseparable state pb s - This 
be fixed by introducing local ancillas on each 
as labels that indicate which term 
defining p hs 



can 

system acting 

in the decomposition is considered 

E t E,c t (K fc tX^tl)i ® •■■ ® (\t,ktXt,k t \) n ® pY® P k t t, 



M ai \xi 
P(a\x) 
correlations 



M 



kt 



tr [M ai \ Xl ® . . . <g> M an \ Xn p bs ] . 



of the form (A4) can be 



Eq. (A2| 



we finally get 

and thus any 
written as in 



Let Qt denote the set of quantum correlations defined 
for the parties belonging to t and let Qb B denote the set 
of biseparable quantum correlations, which is a subset 
of the n-partite quantum correlations. From (A3), it is 
clear that Q bs is convex and that its extreme points are 
of the form P^ xt (a t \x t )PQ Xt (at' \x t >) where PQ Xt (a t \x t ) is 
an extremal point of the set Q t . 



Appendix B: DIEWs for genuine n-partite 
entanglement 

Here we derive the maximum value that the expres- 
sion J n , see Eq. (5), can take for measurements made 
on a biseparable quantum state (the biseparable bound) , 
and for Svetlichny models (the Svetlichny bound). From 
(A3), it follows that the biseparable bound is obtained 



by maximizing /„ over all correlations of the form 

P (*\x) =EE P Q ( 5 *N P Q I s *') ■ ( B1 ) 

In a Svetlichny model, on the other hand, one considers 
correlations of the form 



P(a\x)=J2J2 Pkt (*t\xt)P kt fa \x* ) 



(B2) 



i.e. no constraints (apart from positivity and normal- 
ization) are imposed on the joint terms P kt (at\xt). But 
more refined models a la Svetlichny can be introduced, 
see e.g. |21j . the more constraining one being the one 
where the joint terms P kt {at \xt) are assumed to be no- 
signalling. We will compute the bound on I n assuming 
these no-signalling constraints and see that even in this 
case there exists a gap between the quantum biseparable 
and Svetlichny bounds (note also that the no-signalling 
bound that we will derive actually coincides with the 
more general unconstrained Svetlichny bound because I n 
is an inequality involving only full n-partite correlators) . 

We will now prove the biseparable bound given in 
Eq. ^ for n > 3 using induction on the number of parties 
n. Let us start with the first step of the induction, which 
consists of showing that the inequality holds for n = 3. 
By linearity of /3 in the probabilities, convexity of the de- 
composition Q for tripartite biseparable quantum cor- 
relations and the fact that J3 is invariant under any per- 
mutation of the parties, it is sufficient to prove the bound 
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for correlations of the form PQ(ab\xy)PQ(c\z). Moreover, 
it is sufficient to consider the case where Pq(c\z) is ex- 
tremal, i.e., one where every c is determined unambigu- 
ously as a function of z. To this end, let us label the eight 
distinct deterministic assignments of outcome c for given 
input z by 7 = (71, 72, 73) € {±1} 3 , where the determin- 
istic (quantum) probability distributions are such that 
P Q (c\z) = P 7 (c\z) = 1 if 7 Z = c and P 7 (c|z) = oth- 
erwise. Substituting these eight possible strategies into 
I3 [see ([5])] gives eight different bipartite Bell expressions 
for the system AB: 
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k=2 



(B3) 



Efe=n+i /fc-n-E^+i, we obtain 



In+l ~ 



E E - 

£C„+1 = 1 a„ + i= — 1 



-ri+i-PQ(a„+i|a; n+ i)/^" +1 (a n+ i,x„ + i) 



E3n+3 r pk-j 
k=n+l Jk-n-^n 



(B5) 



where we have defined I n 

EtnA-n+j^n (note that E k is different from zero 
only if n < k < 3n) and write I n (a n+ i, x n +\) 
to remind that I n is evaluated on the distribu- 
tion PQ(ai i ... ! k\xx ! ... < k)PQ(dk+i,..., n \xk+i,..., n ,a n+ i,x n+1 ) 
which depends on a n+ i and x n +i- Using the fact that 
a n+ i — ±1 and that the probabilities P(a n +i|x„+i) are 
bounded between and 1, it follows from (B5l that 



where g 7 (k) — J2 z =i Izfk+z-z- Determining the bisepa- 
rable bound of I3 therefore amounts to finding the max- 
imum of the above Bell expressions over arbitrary bipar- 
tite quantum correlations Pq{ab\xy) = tr[M | x ® M b \ y p], 
i.e., finding the corresponding Tsirelson bounds. It can 
be easily verified that the eight expressions are either 
vanishing or equivalent, under relabelling of inputs and 
outputs, to two times the 3-input chained Bell inequal- 
ity [22]. The biseparable bound on I3 is thus equal to 
the two times the Tsirelson bound of the 3-input chained 
Bell inequality [23], i.e., I3 < 6v3, which is in agreement 
with Eq. (pji 

Next, we need to show that whenever inequality ([5| 
holds for n parties, it must also hold for (n + 1) par- 
ties. By linearity of /„ in the probabilities and convex- 
ity of the decomposition (BI), it is sufficient to com- 



pute the bound over all product correlations of the 
form PQ{at\xt)PQ{a,f\xt) . Moreover, since I n +i is sym- 
metric under any permutation of parties, it is suffi- 
cient to consider the biseparations t — {1, 2, . . . , k}, 
t' = {k + I, . . . , n + 1} where I < k < [f J • With 
obvious notation, we write the corresponding corre- 
lations as PQ(ai i ... )fc |a;i ) ... >fc )PQ(a A+li ... i7l+1 |x fc+ i ) ... iT , + i). 
Bayes' rule and the fact that quantum correlations sat- 
isfy the no-signaling principle allow us to write the 
quantum correlations for parties k + I to n + 1 as 

PQ( a k+l,...,n\%k+l,...,n+li a n+l )Pq (<*n+l pfc+1,.. . ,n+l ) = 

PQ{ak+i,..., n \xk+i,...,n,a n+ i,x n+ i)PQ{a n+ i\x n+ i). With 
straightforward algebra, it can be seen that the quantity 



71 +1 



rpk 

E n +1 

bution can be rewritten as 



k)E{x) for this type of distri- 
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E n+i = E E a n+ iP Q {a„ + i\x n+ i)x 

x„+i=la„+i=-l ( B4 ) 
E k n - x ^{a n+1 ,x n+1 ) 

where the n-partite quantity En Xn+1 (a n+ i, x n +i) 
depends on the values of a n +i and x n + 1 
since it is evaluated on the distribution 
PQ(ai,...,k\ x i,.~,fs)PQ{ a k+i,... 

Inserting this expression in the definition I n +i = 



In+l < E l / n" +1 ( a ™+l,a;n+l) 



(B6) 



Let max \I n \ — max i/^ denote the maxi- 
mal value of the quantity ±7^ taken by any 
biseparable correlations. Since the distribution 

PQ(ai,...M x h---,k)PQ( a k+l,...,n\Xk+l,...,n,an+lX n +l) IS 

biseparable, we clearly have that \I n n ~ 1 (a n +i, x n +i)| < 
max H n , independently of the values of a n +i and x n+ \. 
We thus find 



-1 < E max±/ « 



(B7) 



An important point to note now is that ±I n for all 
j = 1,2,3 are equivalent to /„, i.e., one can obtain the 
expression for /„ by starting from any of ±I n and ap- 
plying a different labeling for the inputs and/or outputs. 
To see that this is the case, we first note from the def- 
inition of fk that — lf % = I n . In other words, we can 
obtain 1^ by starting from I n and applying the mapping 
a n — > —a n to each of the Clearly, the same argument 
also demonstrates the equivalence between I n and —I n . 
Next, note that if the following mappings are applied to 
the definition of E k , namely X fi r X y, = (x n mod 3) + I 
and a n — > a! n — a n except for x n — 3 in which case 
a' n = — a n , then we can also obtain the expression for I n 
from I n +1 , thus showing their equivalence. We thus have 
that maxifi < max/„ and thus 



In+i < 3max/ n 



(B8) 



Now, by the induction hypothesis we know that max /„ < 
2 x 3"~ 3 / 2 for all biseparable correlations . It then follows 
from Eq. |B8]) that I n+1 < 2 x 3™" 1/2 , which completes 
the proof. 

Note that, following the above reasoning, any witness 
for n parties and m inputs per party which can be writ- 
ten in the form J^k fkEn> where fk '■ Z — > K is a function 



satisfying f k+ r 
ties. 



±/fe, can be generalized to more par- 
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Note also that the above biseparable bound for the 
witness ^ is tight, as it can always be achieved by per- 
forming the following local measurements 

{ cos (j>(xi)a x + sin (f>(xi)a y : i — l,...,n— 1 
g z : i — n 

(B9) 

with 4>(xi) = (^3^ + g( n 1 _i) ) 7r on the biseparable state 

\fp n ) = \GHZ„^i)»\0). (BIO) 

Finally, let us note that the same procedure as the 
one detailed above can be followed to obtain the (no- 
signalling) Svetlichny bound of I n +\- The only differ- 
ence is that in the first step of the induction proof, we 
must compute the no-signalling bound of the J3 inequal- 
ity instead of the quantum biseparable bound. As in 
the quantum biseparable case, this reduces to computing 
twice the (bipartite) no-signalling bound of the 3-input 
chained Bell inequality, which gives 7 3 <2x6 = 12 and 
thus /„<4x 3™ -2 . To show that this bound is tight con- 
sider a strategy where the n-th party is separated from 
the rest and always outputs 1. Using an analysis similar 



to the one leading to Eq. ( B6 1 , we can then see that the 
remaining (n — 1) parties are playing an effective game 
defined by I^-i + In-i + ^n-i which can be shown to be 
equivalent to twice the expression J n _i. It is a simple 
exercise to show that the algebraic maximum of I n —\ is 
2 x 3"~ 2 and that this is always achievable by (n— 1) play- 
ers that are constrained only by the no-signalling princi- 
ple. Therefore, with this particular strategy, one achieves 
4 x 3™~ 2 , which saturates the bound derived above. 



Appendix C: Biseparable model for projective 
measurements on the tripartite GHZ state 

Here we present a biseparable model that simulates von 
Neumann measurements on noisy tripartite GHZ states 

p = y loro + iii X Q(X) + iii| + i (C1) 
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of visibility V = \ . Clearly this allows to simulate states 
with V < \ as well, by mixing this model with another 
one in which all parties produce uniformly random out- 
comes. 

The model is presented as a protocol in which a source 
distributes quantum states and random variables to the 
parties. All parties can share the random variables, but 
since only biseparable states may be used in the model, 
not more than two parties at a time can share a quan- 
tum state. After distribution, the parties receive their 
respective measurement directions x, y or z belonging 
to the Bloch sphere, and measure their quantum system 
and/or process the information they received accordingly 
in order to produce binary outcomes A, B, C = ±1. 

We check that the outcomes produced by the mo del are 
identical to the ones found when measuring state (CI I in 
the given bases. 



1. The model 

Before the parties receive their inputs, a common 
source chooses a party p at random, say Charlie (p = C), 
and sends him the vector 



A = (sin a cos /?, sin a sin /3, cos a), 



(C2) 



uniformly chosen on the sphere S 2 . The source then also 
provides the two other parties, Alice and Bob in this case, 
with the quantum state 



\<f> AB ) =cos-|00)+sin|e-^|ll) 



(C3) 



Moreover, the source sends to all parties the signs 
sab, sac, s bc = ilj which are independently and iden- 
tically distributed with Prob(s = +1) = ^±^§. 

At the time of measurement, Alice and Bob measure 
their system according to x and y and get outcomes a,b — 
±1, while Charlie calculates c = sg(A • z). The parties 
then output respectively A = s A B s Ac a , B — sabSbc^ 
and C = sbcsacc 



2. Correlations obtained by the model 

Here we compare the correlations that are created 
when Alice, Bob and Charlie apply the biseparable model 
above to the ones that they would get by measuring 
state (CI ). Since the model treats equally each party, its 



correlations are symmetric under exchange of the parties 
and we only need to check the following three relations: 



(A) = 
(AB) = i cos 9 X cos 9 y 



(C4) 
(C5) 



(ABC) = — smQ x s'md y sm6 z cos((p x + <p y + ip z ). (C6) 

Here we parametrized the parties' measurements in terms 
of spherical coordinates: 



x = (sin 6 X cos (p x , sin 9 X sin tp x , cos ( 



(C7) 



and similarly for y and z. 

For definiteness, in the following we use the brackets 
(q) to express the expectation value of a quantity q with 
respect to a quantum state, and the bar 



1 

■in 



2tt 



df3 / da sinaf (a, /3) 



(C8) 



for the average of a function f(a,/3) over the random 
variable a and j3. We always start by considering that 
Charlie is the special party chosen by the source, to 
whom the hidden vector A is sent, and perform the 
symetrization afterwards. Symmetrized quantities q^ — 
|(<7 P= A + <Zp=B + <Zp=c) are indexed by the symbol 0- 

For simplicity we average over the depolarization signs 
sab, sac, s bc only at the end of the calculation. 



a. Single-party expectation value One can easily ver- 
ify that the outcomes of each party is locally random ac- 
cording to the model, in agreement with equation (IC4I). 



b. Two-party correlators Let us first calculate cor- 
relations (ab), (ac) and (be) for the case in which Charlie 
receives the vector A, and average over the choice of the 
party p being alone afterwards. 



The first term is easily found from equation ( C3 1 , 
which gives 

(ab) = cos 9 X cos 9 y + sin a sin 9 X sin 9 y cos(/3 + ip x + tp y ) 

' (C9) 

and thus (ab) = cos 9 X cos 9 y . 

To compute the correlation (ac), it is useful to write 
Alice's state as 



PA = tT B {\$ AB X$ AB \) = cos 2 — |0X0| 



sin 2 f |1X1|- 



(CIO) 

The expectation value for Alice's outcome is then (a) = 
cos a cos 9 X . 

Concerning Charlie, his outcome c is totally deter- 
mined by A and z. For simplicity we assume that 
< 9 Z < 7r/2, but the other situation can be treated 
similarly. In this case, for a < f — 9 Z , Charles al- 
ways has c = +1, and for a > 5 + 9 Z) he always has 
c = — 1. Now for a £ [~ — 9 Z , | + 9 Z ], we can write 
X_l_ = [— $ c + ip z , $ c + cp z ] the interval of f3 for which the 
product A • z > 0, and X_ = [ip z — Tr,ip z + tt]\T + its com- 
plement, with $ c — arccos(— cot a cot We thus have 
the three following cases: 



1. a < 



9 Z . In this case c 



(ac) — cos a cos ( 



2. a > 

(ac) -- 



d. 2 



■f 9 Z . In this case c 
cos a cos 0, . 



-1, which gives 
T, which gives 



<a<f 



In this case one has (ac) 



cos a cos 6» x (xi + (/?) - Xi- (P)) where 
XrC8) 



1 if P £ 1 
if P £ 1 



(Cll) 



is the indicator function. 



In total, after integration over a and ft, this gives 

(ac) = \ cos cos 9 Z . 

Similarly one can check that (be) — \ cos 9 y cos 9 Z , and 
so once averaged over the choice of the party being alone, 
the bipartite correlations are given by 



(ab) 



= — cos 9 X cos ( 



(C12) 



Finally, we need to apply the signs sab, s ac> s bc- 
Since the expectation value of the product of two signs 



is (sab sac) — j > it follows that overall the correlation 
between two outcomes produced by the model is 



(AB) = - cos# x cos 6^, 



(C13) 



in agreement with equation ( C5 ) 



c. Three-party correlators We now proceed to cal- 
culate the tripartite correlations that arc created by the 
model. Let us consider the three preceding cases sepa- 
rately again: 



1. a< | 



9 Z . In this case c = +1, so (abc) — (ab) as 



given by equation ( C9 ) 



2. «>f 



9 Z . In this case c : 



T, so (abc) 



(ab). 



3. f — 9 Z < a < | + 9 Z . In this case one has (abc) = 
(ab)(xx + (P)-Xx.(P)). 

In total, the tripartite correlation are found after inte- 
gration over a and p to be 



(abc) = 



1 



sin o ~ sin ( 



: cos((p x + ip y + ip z ), (C14) 



which does not depend on which party p is alone. 

One can check that the application of the signs 
has no influence on the tripartite correlations since they 
cancel out, and so (ABC) = (abc). The model thus re- 
produces the expected correlations (C6). 



Appendix D: Characterizing biseparable correlations 
through SDP in the tripartite case (n = 3) 

It has been shown in IS] how to define a hierarchy 
of SDP that characterizes the set Q2 of bipartite quan- 
tum correlations. Note that biseparable quantum cor- 
relations P(abc\xyz) can be written as a finite sum of 
such bipartite quantum correlations using the fact that 
each P k (c\z) in Q can be taken to be extremal and us- 
ing the fact that there are a finite number of such ex- 
tremal points corresponding to classical, deterministic 
strategies. It therefore follows that biseparable quantum 
correlations can be characterized using a finite number 
of the SDP hierarchies introduced in [TH]. Note however 
that the number of single-party deterministic strategies, 
and thus the number of terms in Eq. Q, grows exponen- 
tially with the number of measurement settings, making 
such an approach impractical even for small problems. 
Here we introduce an alternative SDP approach that has 
better scaling properties. 

Our approach is based on the observation that the ten- 
sor product separation pab ® Pc a t the level of states 
can be replaced by a commutation relation at the level of 
operators. Specifically, let s = {AB/C,AC/B,BC/A} 
denote the three possible partitions of the parties into 
two groups. Then, P(abc\xyz) are biseparable quan- 
tum correlations if and only if there exist three arbitrary 
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(not necessarily biseparable) states p s and three sets of 
measurement operators {M^ 
holds, i.e. 



Ax ,M% ly ,M° lz } such that f6 



P{abc\xyz) = ]T tr[M Q % M b s |y ® M c s |z p s ] , (Dla) 

s 

where measurement operators corresponding to an iso- 
lated party commute, i.e., 



K, c/S '<?' /B ] = ° 

BC/A 



(Dlb) 



r M ^, A M BC/A 



]=o. 



Here, for 



convenience, the normalization factor of the 
has been absorbed in the quantum states 



RHS of ( [Dla 
/''• 

Let us start by showing that any correlations of the 
form Q is of the form (Dl). Note first the well- 



known fact that if the dimension of the Hilbert space 
is not fixed, any single-party correlations P k (c\z) admits 
a quantum representation P k (c\z) = tr[M^ z pp] where 

the operators Mh x commute among themselves. Indeed, 

let c = (ci,...,c m ) and P(c) = YiT=i P( c z\z), then a 
quantum representation of P k (c\z) is achieved by defin- 
ing p h c = £ E P(c)|c)(c| and m\ z = £ s <5(c z , c)\c)(c\. We 
can thus write in Eq. (jij), e.g., J2k PQ(ab\xy)P k (c\z) = 

tv[M a \ x ®M b \ v ®M c \ z (pAB®Pc)} with the operators 
commuting between themselves. The decomposition p} 



is thus clearly a particular case of (Dl ). 

Let us now show the converse. Let us suppose first 
that s — AB/C and let c = (ci,...,c m ). Since 
the projectors M*, z commute, the projectors M(c) = 

Y\T=i ^c z \z defines a valid measurement and thus pc = 
trc [M|p s M|] is a proper (unormalized) state of sys- 
tem AB. Define P^(ab\xy) = tr[M* 



2 | X 



M b\yPc\ and 



P c {c\z) as the deterministic point satisfying P c (c\z) = 
<5(c z ,c). It is then easy to see using the properties of 
the operators M° ]z that ti[M^ x ® M^ y ® Mf |z p s ] = 



Ee*(p«ic)tr[M | x ®Af 6 | 1 ,®Mc/!> a ] = £ s <5(c z , clt^M^® 
-^b|yPc] = PQ{ a b\ x y)P c { c \ z )- This last expression is 
of the same form as the first series of terms in Q (since 
any deterministic point P c {c\z) can be seen as a single- 
party quantum point Pq{c\z)). A similar argument for 
the other values of s implies that any correlations of the 
form (Dl ) are of the form Q. 



We thus have reduced the problem of determining 
whether given correlations P(abc\xyz) are biseparable to 
the problem of finding a set of operators satisfying a finite 
number of algebraic relations (the projection defining re- 
lations of the type M' |s M« |a = 8 a , a ,M° , £ a M* = 1 



and the commutation relations (Dlb)) such that (Dla 
holds. Such a problem is a typical instance of the SDP 
approach introduced in 19J, which generalizes the results 
of HHJ. Specifically, it follows from the decomposition 



(Dl) and the results of 19J that it is possible to define 
an infinite hierarchy of sets {Q2/1 ■ 3 ' ~ ^ ■ • ■ ' °°} with 
the following properties: i) deciding whether given cor- 
relations belongs to Q 3 2 /i can be determined using SDP 
(though the size of the SDP increases with j); ii) the 
sets Q J 2 /i approximate better and better Q2/1 from the 
outside, i.e., Q2/1 C 



H,j and Qi,^ ^= Qyi for au 3- If 



given correlations are found not to belong to some Q2/1' 
we can thus conclude that they do not belong to Q2/1 
and that they reveal genuinely tripartite entanglement. 
Further, it is possible in this case to derive an associated 
DIEW from the solution of the dual SDP. Note that mod- 
ulo the assumption that the tensor product between the 



three measurement operators in (Dla) is equivalent to the 
commutation of these operators — an assumption that is 
true in finite-dimensional Hilbert space, but which still 
represents an open question in an infinite-dimensional 
Hilbert space — it can be shown that the hierarchy of 
relaxations converges to the set of biseparable states, i.e., 
linij-j.oo Q2/1 = *?2/i- Thus if given correlations are not 
biseparable this will necessarily show-up at some finite 
step in the hierarchy. 



